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f~ i . Abstract. Different approaches to the problem of a gauge- invariant operator definition of 

transverse-momentum dependent parton densities (TMDs) are reviewed and compared from the 
point of view of their compatibility with the operator definition of coUinear (integrated) parton 
^j^' distribution functions. In particular, geometry of the longitudinal and transverse gauge links 

' (Wilson lines) both in unsubtracted quark TMDs and in the associate soft factors is addressed. 

A possible connection between different operator definitions is also discussed. 



> , 

. Classical deep inelastic IH — t- I'X scattering experiments deepened considerably our 

I knowledge of nucleon structure. Measuring the momentum of the outgoing lepton /' we learn 

■ about the distribution of partons inside the nucleon. In a reference frame where the nucleon 

, moves (infinitely) fast, this information is accumulated in the parton distribution functions 

(PDFs) faix^Q"^) in terms of the partonic degrees of freedom: e.g., the Bjorken variable xbj 
relates to the fraction of the longitudinal momentum P of the parent hadron possessed by a 
parton of the flavor a. Such collinear (integrated) PDFs can be properly defined as completely 
gauge invariant (nonperturbative) hadronic matrix elements (here and in what follows, the light- 
^ . cone components of four- vectors are = {p^ ibp^)/-v/2) 



fa{x,fl^) =27^ ^ (P \Mr,0^)[r,0-]n7^MO-,0±)\ p) (1) 

with renormalization-group properties controlled by the DGLAP evolution equations (for review 
and Refs. see [1]). The Wilson line ['^ , 0^]„ will be introduced below. Moreover, one can 
relate the moments of the collinear PDFs = jdx x^~^ fa{x) to the matrix elements of the 
local twist-two operators = {p~^)~^ {p \^-il^a{0){'y^iD^ ...iD~^}sym.'ipa{0)\ p) arising in the 
operator product expansion on the light-cone, thus making them well-defined objects from the 
field-theoretical point of view . 

Another attractive feature of the collinear PDFs ([1]) is that they imply a clear interpretation 
as the probability for a parton inside the nucleon to have the longitudinal momentum /ciong. = 
xP\ong.- This interpretation is most naturally established when QCD is canonically quantized 
(and subsequently renormalized) on equal- "light-cone-time" surfaces ^'^ = in a class of singular 
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non-covariant gauges [3]. This parton number interpretation holds, in the hght-cone gauge 
= 0, in higher order calculations as well [UEllG]. However, an important problem must be 
solved in order to make the calculations in the light-cone gauge feasible. Even after imposing 
the gauge condition = 0, the gauge in not completely fixed: one may still perform a 
independent gauge transformation U{^^,^_l) by virtue of d^U{^,^ ,^±) = d / dx~ U {^~^ , = 0. 
Therefore, fixing the gauge is equivalent to imposing certain boundary conditions on the gauge 
field (see, e.g., [SI El El E] ) . Going over to the conjugate momentum space, one observes that the 
ambiguity in the behavior of the gauge field at light-cone infinity — )• oo maps over ambiguity 
of the gluon Green function at small — )• 0. A key issue is, therefore, how to get rid of extra 
complications due to the emergent "spurious" singularity ~ [^7"*"]"^ in the free gluon propagator 

The uncertainty of the pole prescription in Eq. ([2]) corresponds to the residual gauge freedom 
and can be treated without changing the gauge-fixing constraint = 0. 

There are several possible pole-prescription-fixing procedures that are compatible with the 
light-cone gauge and have been shown to give correct results (at least, up to the O(a^)-order). 

In particular, the principal value prescription ^^^^v = lini,j->.o ^ + g+^-i-q ) used 

in Ref. [lOj to evaluate the DGLAP kernel in the next-to-leading order. Non-symmetrical 
advanced and retarded pole prescriptions are also possible [3 [8]. Although these methods work 
in some situations, the only pole prescription which is consistent with the equal-time canonical 
quantization in the light-cone gauge is the Mandelstam-Leibbrandt one [llj : 

' lim — — j = lim — -^^ 3- — = lim — — — - , (3) 



[q+] v^o [q'^]ML v^o {q ■ n+){q ■ n ) + ii] v^o {q-n ) + ir]{q ■ n+) 

where and in what follows are the light-like vectors (n^)^ = , n+n^ = 1, and = 
means equality in the sense of the theory of distributions. It was shown that the free gluon 
propagator supplied with the ML pole prescription can be directly derived following the equal- 
time quantization procedure and is compatible with well-established results at least up to 0{as)- 
order [ll|8l[T2]. The main difference between the ^"-independent prescriptions and the ML one 
([3]) originate in the different situation of poles in the q^ plane, as it is shown and explained in 
Fig. 1. Thanks to this feature, one can perform a Wick rotation of the integration contour to 
the Euclidean momentum space, and the ultraviolet divergences can be analyzed by means of 
the usual power counting procedure in the Euclidean space. This observation anticipates the 
absence of overlapping divergences in the loop calculations with the ML prescription. 

In is worth noting that albeit the pole-prescription issues mentioned above may reveal 
themselves in the course of the calculation, they are non-visible in the case of the collinear PDFs, 
by virtue of the cancelation of the soft divergences in the virtual and the real gluon exchange 
graphs. However, those issues are crucial and unavoidable in unintegrated PDFs [2j, which are 
introduced in the factorization approach to the semi- inclusive processes. The picture of the 
nucleon revealed in the DIS experiments, being essentially one-dimensional, is still incomplete: 
in this "longitudinal" picture the transverse degrees of freedom of the partons are eliminated by 
definition and the 3D-structure remains inaccessible. The study of semi-inclusive processes, such 
as semi-inclusive deep inelastic scattering (SIDIS), the Drell-Yan (DY) process, hadron-hadron 
collisions, or lepton-lepton annihilation to hadrons, where (at least) one more final or initial 
hadron is detected and its transverse momentum (and, possibly, its spin) is observed, calls for 
the introduction of more involved quantities — unintegrated transverse-momentum dependent 
(TMD) distribution and fragmentation functions (see [T^ and Refs. therein). Moreover, a 



Wick rotation 



Figure 1. Location of the poles in the complex q plane: the poles of the light-cone gluon 
propagator with the ML prescription (1) and the poles of the propagator in a covariant gauge 
(2) belong to the same quadrants, so that the clock-wise Wick rotation is allowed. The poles 
of the light-cone propagator with the principal- value prescription (3), in contrast, impede that 
rotation. 



variety of applications of the TMDs approach has been found in the phenomenology of polarized 
hadronic processes. It was conjectured and corroborated that the idea of TMD parton densities 
inside the nucleons is directly applicable in the theory of single-spin asymmetries (see, e.g., 
|14j and Refs. therein). Recently, interesting possibility that linearly polarized gluons inside 
unpolarized protons can affect the cross-sections of the scalar and pseudoscalar Higgs boson 
production in the gluon-gluon fusion channel was proposed and discussed within the TMD 
approach in the Refs. [T5] . 

In the present paper, I give a brief account of how the problem of the emergent singularities 
beyond the tree-approximation is approached in different operator definitions of the (quark) 
TMDs. At the one-loop level, the following three classes of singularities are expected: {i) simple 
ultraviolet poles which must be removed by the standard renormalization procedure; {ii) pure 
rapidity divergences, which depend on an additional rapidity parameter, but do not jeopardize 
the renormalizability of the TMDs, and can be safely resumed by means of the Collins-Soper 
equation; {Hi) highly undesirable overlapping divergences: they contain the UV and rapidity 
poles simultaneously and thus break down the standard renormalizability of TMDs, calling 
for a generalized renormalization procedure in order to enable the construction of a consistent 
operator definition of the TMDs. Before getting started with the analysis of the divergences, let 
us try to learn something from the tree-approximation. The simplest "unsubtracted" definition 
of "a quark in a quark" TMD ([A]n means that we use the light-like longitudinal Wilson lines), 
which allows a parton number interpretation in the light-cone gauge, reads 

x[oo",^±;oo",oo±]l7^[oo",oo_L;oo~,0_L]i[oo",0^;0",0_L]nV'a(0",0_L)| (4) 

with ^+ = 0. Generic semi-infinite path-ordered gauge links evaluated along a given four-vector 
w are defined as [oo;.^]^^. = "Pexp [—ig f^dr Aat°'{^ + wt)\ , where, in the case under 
consideration, the vector w can be either light-like wl = , (n^)^ = 0, or transverse wt = I- 



Although the formal integration of definition ^ over k± yields the collinear PDF, Eq. ([T]), 
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{P \Mr,0^)[r,0-]n7+MO~,0±)\ p) = fa{x) , 

(5) 

this is only well-justified in the tree approximation, because the rapidity divergences in the loop 
corrections prevent a straightforward reduction. I would like to emphasize that the normalization 
of the above TMD 
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can be obtained following the canonical quantization procedure in the light-cone gauge, where 
the longitudinal Wilson lines disappear and the equal-time commutation relations for the quark 
creation and annihilation operators {a^(A;, A), a(fc. A)} lead immediately to the parton number 
interpretation of the TMD: 
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Use of the "tilted" gauge links doesn't obey this requirement. 
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Figure 2. Geometrical structure of integration contours in the unsubtracted TMDs with 
the light-like (upper panel) and off-the-light-cone (lower panel) longitudinal gauge links and 
their symbolic reduction to the collinear PDFs. In the former case, the transverse gauge links 
cancel completely after the fc^-integration, while the longitudinal gauge links reduce to the one- 
dimensional light-like connector [^",0"]. In the off-the-light-cone situation, the cancelation of 
the transverse gauge links at infinity is not, at least, straightforward. Moreover, the integrated 
configuration contains two non- vanishing off-the-light-cone gauge links which are not equivalent 
to the simple connector [^^,0^]. Beyond the tree-level, the renormalization group properties 
of those two objects are also different. I put the interrogation marks next to the transverse 
gauge links at infinity since I'm not aware of any consistent treatment of them in the TMD 
formulations with off-the-light-cone (tilted) Wilson lines. In contrast, the transverse gauge links 
appear naturally in the "light-cone" frameworks. 




Figure 3. Comparative geometry of the Wilson lines in unsubtr acted soft factors and 
visualization of the reduction to the collinear case. Upper panel shows the soft factor in the 
momentum space, as proposed in the Refs. [19\ I21j . Lower panel presents the tilted off-the- 
light-cone integration paths in the impact parameter space, as well as the result of the reduction 
to the collinear b± ^ configuration. 




Figure 4. Comparative geometry of the Wilson lines in the subtracted soft factors. Upper panel 
corresponds to the soft factor of the TMD distribution function which enters the factorization 
formula ([9]). Lower panel shows the longitudinal gauge links shifted off the light-cone, which are 
used in the factorization approach (jlOp . 

Going beyond the tree-approximation, one encounters a bunch of singularities mentioned 
above. To overcome the problems related with them, different frameworks have been proposed. 

Adopting the covariant Feynman gauge, Ji, Ma and Yuan developed a framework which makes 
use of the tilted (off-the-light-cone) longitudinal gauge links lined up along the vector / [16]. 
In a covariant gauge, the transverse gauge links at light-cone infinity cancel, and the rapidity 



cutoff C = (2p • 'ri'B)'^ /\'n\\ is introduced to control the deviation of the longitudinal gauge links 
from the light-like direction. A subtracted soft factor contains the non-light-like gauge links as 
well. Within this approach, the off-the-light-cone unsubtracted TMDs where the light-like vector 
n~ in the longitudinal Wilson lines is replaced by the tilted vector ne = (— e^^^, l,Oj_) do not 
satisfy the relation ([5]), even in the tree approximation — cf. Fig. 2 and the caption. However, 
one can design a "secondary factorization" method which allows the expression of this TMD 
(transformed to the impact parameter space in terms of a convolution of collinear 

PDFs and perturbative coefficient functions at small b_L [TB]. Another subtraction method, also 
in covariant gauges, but without explicit off-the-light-cone regularization in the unsubtracted 
TMD, was developed in Refs. |17[ I18j. The corresponding geometry of the light-like and tilted 
Wilson lines in the soft factors is shown in Fig. 3, lower panel. 

In our works [TU] we proposed to study the renormalization-group properties of the 
unsubtracted quark TMD (jj]) and to make use of its one-loop anomalous dimensions in order 
to reveal the simplest minimal geometry of the gauge links in the soft factor which allows one 
to get rid of the mixed rapidity-dependent terms. We showed (in the leading O(as)-order) 
that the extra contribution to the anomalous dimension is nothing but the cusp anomalous 
dimension [20l [12] . Note that in these works we adopted the light-cone gauge supplied with the 
^"-independent pole prescriptions. In subsequent works we showed that our approach works in 
the case of the ML pole prescription as well, Refs. [21], and that it can be consistently used to 
formulate a generalized definition of the quark TMD with a non-minimal spin- dependent term 
in the Wilson lines, Refs. p2] (see also [23]). 

Compared to the Ji-Ma-Yuan approach, we followed a different strategy. Making the 
assumption that the parton number interpretation ([7]) must hold for TMDs in the light- 
cone gauge (like it holds in the collinear PDFs), we are in a position to derive a gauge- 
invariant operator definition of the TMD. In other words, starting from the requirement of 
the probability interpretation in the light-cone gauge and adding, step by step, the minimally 
necessary gauge links, we would end up with a gauge invariant operator definition of the TMD 
without undesirable singularities [23]. The generalized definition of the quark TMD we proposed, 
reads 

^ SF{n+,n--e)L-p^{n+) ^ ' 

where the soft factor Sp and the self-energy factor Lp are defined in Refs. [191 [26] , see also 
Figs. 3 and 4. The rapidity regulator is ^ = ^"^^ . 

Let me suggest the following conjecture concerning the generic structure of divergences in 
(any reasonable) operator definition of the TMD beyond the tree-approximation: 

The contribution of the overlapping ~ l/e®\\i9 singularities to the renormalized TMD can 
be expressed either in terms of a finite number of the cusp anomalous dimensions which are 
known in the theory of Wilson lines/loops up to the 0(ag)-order — in this case, their treatment 
consists of the subtraction of (a finite number of) corresponding cusped soft factors] or those 
singularities depend on the degenerate rapidities ~ In^o = In(n^)^ — in that case, one has to 
subtract the self-energy soft factors which consist, in contrast, of the "smooth" infinite gauge 
links without any obstructions (cusp or intersections). The conjecture is, therefore, that there 
is no other sort of unphysical singularities in the loop corrections to the TMDs in any order of 
as- 

In the leading order, we demonstrated the validity of the above statement in our works 
\19 \ \21\ [25l 126] . Recently, some interesting aspects of the TMD factorization, light-like Wilson 
lines and renormalization properties of the PDFs have been studied adopting the methodology 
of soft-collinear effective theory (SCET) in Refs. [27]. An important question remains, however, 
how our TMD can be built into an appropriate factorization formula for semi-inclusive hadronic 



tensor. In our approach, the following factorization formula is supposed to be valid 



' ' SF{n+,n-;e)L-\n+) SD{n+ ,n--e)L-\n") ' ^ ' 

where the geometrical structure of the soft factor is consistent with that of the collinear PDF 
\12\ [28] and doesn't break the number interpretation. The explicit proof of the conjectured 
factorization and of absence of double-counting is in progress. 

Recently, Collins proposed a new definition of the (quark) TMD [29] which is built into the 
factorization formula for the semi-inclusive hadronic tensor (up to power corrections) 

W'^'' = \H{Q, f,fr ■ (x, kr,f^, Cf) ® {z, zk'^;f,, (d) + ■•. , (10) 

where all soft factors are absorbed into the TMD distribution and the fragmentation 

functions, so that there are no separate soft factors in factorized structure functions, e.g.. 



Here the soft factors depend on the light-like or the tilted nA,B vectors (for details, 
see [29]). Note that the TMD ([TT]) is defined in the impact parameter space T{x,b±) = 
JcPk± Q-^^i-^i- T{x,k^), so that it is, in fact, a "semi-integrated PDF" and the reduction 
to the collinear case corresponds to the limit b± — )• 0. Some phenomenological and lattice 
applications of this approach have already been discussed in Ref. [30]. The geometry of the 
gauge links in the soft factors is presented and explained in Fig. 4. 

Several open questions are still to be answered: (i) How should one prove the complete 
gauge invariance of the TMD (|lip ? It is formulated in the covariant Feynman gauge where the 
transverse gauge links at light-cone infinity vanish. What will change, if we adopt some physical 
(axial) gauge? (ii) In particular, how one should treat the T— odd effects in the axial gauges 
given that the structure of the transverse gauge links at light-cone infinity is not yet clarified 
in the TMD ([U])? (in) After reduction to the cohinear PDF (in the case of TMD ([n]), this 
corresponds to the limit b± ^ 0), there is neither a mutual compensation of the longitudinal, 
nor that of the transverse gauge links (if introduced in the usual manner), see Fig. 4. Hence, 
the geometrical structure of the gauge links in the collinear PDF obtained from the TMD (jlip 
seems too cumbersome to be simply included in the standard DIS factorization scheme (see, 
e.g., Ref. [I2l[28]). 

To conclude, several approaches to the problem of the factorization of semi-inclusive processes 
which make use of the unintegrated TMD parton distribution functions have been proposed 
and developed so far. There is no a priori clear relationship between these frameworks: e.g., 
corresponding operator definitions of the TMDs may, in principle, describe different objects 
with difi^erent (renormalization- group, gauge invariance, evolution) properties. It is a matter of 
further study to work out these issues. 
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